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ABSTRACT Quantitative trait locus (QTL) mapping studies of Pacific three-spined sticklebacks {Gaster- KEYWORDS 
osteus aculeatus) have uncovered several genomic regions controlling variability in different morphological body shape 
traits, but QTL studies of Atlantic sticklebacks are lacking. We mapped QTL for 40 morphological traits, Gasterosteus 
including body size, body shape, and body armor, in a F 2 full-sib cross between northern European marine aculeatus 
and freshwater three-spined sticklebacks. A total of 52 significant QTL were identified at the 5% genome- genetic linkage 
wide level. One major QTL explaining 74.4% of the total variance in lateral plate number was detected on map 
LG4, whereas several major QTL for centroid size (a proxy for body size), and the lengths of two dorsal morphology 
spines, pelvic spine, and pelvic girdle were mapped on LG21 with the explained variance ranging from QTL 
27.9% to 57.6%. Major QTL for landmark coordinates defining body shape variation also were identified on 
LG21, with each explaining >15% of variance in body shape. Multiple QTL for different traits mapped on 
LG21 overlapped each other, implying pleiotropy and/or tight linkage. Thus, apart from providing confir- 
matory data to support conclusions born out of earlier QTL studies of Pacific sticklebacks, this study also 
describes several novel QTL of both major and smaller effect for ecologically important traits. The finding 
that many major QTL mapped on LG21 suggests that this linkage group might be a hotspot for genetic 
determinants of ecologically important morphological traits in three-spined sticklebacks. 



Understanding the genetic architecture of quantitative traits can pro- 
vide important insights toward elucidating the molecular basis of 
adaptation and evolution (Merila and Sheldon 1999; Erickson et al. 
2004; Phillips 2005; Ellegren and Sheldon 2008). Quantitative trait 
locus (QTL) mapping represents a classical method to uncover geno- 
mic regions controlling variation in phenotypic traits, as well as to 
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gain insights into the distribution of QTL effect sizes of ecologically 
and evolutionarily important traits (Lynch and Walsh 1998; Slate 
2005). Although there are potential problems and biases involved in 
QTL mapping methods {e.g., Beavis 1994; Slate 2005, 2013), they have 
succeeded in revealing surprisingly fine-grained {e.g., Colosimo et al. 
2005) and relatively accurate {e.g., Price 2006) information on the 
location of genetic factors controlling variation in quantitative traits. 

Morphological traits, such as those portraying variation in organ- 
ismal size and shape, typically are complex quantitative traits under the 
control of an interacting network of genes and environmental factors 
(Wu and Lin 2006; Mackay et al. 2009). The genetic basis of morpho- 
logical variation has been studied extensively in domestic animals 
{e.g., Mott et al. 2000; Karlsson and Lindblad-Toh 2008) and plants 
{e.g., Nordborg and Weigel 2008; Bergelson and Roux 2010) but less so 
in wild vertebrates (reviews in: Slate 2005; Slate et al. 2010). One of the 
reasons for this paucity of studies in wild vertebrates relates to the high 
expense and difficulty in raising large numbers of individuals in con- 
trolled conditions. However, such studies have now become possible in 
many emerging model organisms, including the three-spined stickleback, 
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Gasterosteus aculeatus, which can be bred and raised in controlled con- 
ditions with relative ease. 

The three-spined stickleback is emerging as a model system to 
study adaptation and speciation in vertebrates (Bell and Foster 1994; 
McKinnon and Rundle 2002; McKinnon et al. 2004; Gibson 2005; 
Barrett 2010). This small marine teleost fish has repeatedly and in- 
dependently colonized numerous freshwater habitats across the north- 
ern hemisphere since the last ice age and adapted rapidly to new 
environments often in a parallel fashion (Bell and Foster 1994; 
McKinnon and Rundle 2002). The derived freshwater populations 
exhibit morphological phenotypes distinct from their marine ances- 
tors, including changes in body size and shape, armor, and trophic 
morphology (Lavin and McPhail 1986; Bell and Foster 1994; Walker 
and Bell 2000). The genetic bases of these phenotypic traits have 
been the focus of many QTL-mapping studies in the last decade 
(Peichel et al. 2001; Colosimo et al. 2004; Cresko et al 2004; Shapiro 
et al 2004; Kimmel et al 2005; Albert et al 2008; Kitano et al 2009; 
Greenwood et al 2011; Rogers et al 2012). These studies have 
identified several genes and genomic regions contributing to varia- 
tion in quantitative traits and thereby provided important insights 
into the genetic mechanisms of morphological divergence and ad- 
aptation. However, with the exception of two recent case studies, 
which focused on the genetic architecture of a set of correlated traits 
relating to body shape (Albert et al 2008; Rogers et al 2012), most 
of these studies have focused on the genetic architecture of a single 
trait or phenotype (e.g., Colosimo et al. 2004). In addition, almost all 
of the previous investigations have used mapping crosses of Pacific 
origin (i.e., Canada, United States, or Japan; Supporting Informa- 
tion, Table SI), whereas few QTL studies have used populations 
of Atlantic origin (but see Shapiro et al. 2004; Coyle et al. 2007). 
Because QTL can be lineage- and even population- specific (Symonds 
et al. 2005; Chen and Ritland 2013), a comparative approach is 
necessary to distinguish ancestral adaptive variation (i.e., shared 
across lineages) from more recently evolved (i.e., lineage specific) 
adaptations. 

The main aim of this study was to conduct a QTL-scan in three- 
spined sticklebacks originating from the Atlantic lineage, with par- 
ticular focus on ecologically important morphological traits: body size, 
body shape, and armor. To this end, we scanned for QTL associated 
with these traits using 131 microsatellite markers across the whole 
genome in a F 2 full- sib cross generated from marine and freshwater 
populations. We discuss the results in the context of previous research 
done in Pacific three-spined sticklebacks. In addition, we compared 
our genetic linkage map with the physical map generated by BLAST 
searches against the G. aculeatus genome. Apart from adding new 
dimensions to the understanding of the genetic architecture of phe- 
notypic traits in the three-spined stickleback, our findings provide 
insights into the distribution of QTL effect sizes in a number of eco- 
logically important morphological traits. 

MATERIALS AND METHODS 
Sampling and rearing 

One F 2 full-sib family consisting of a total of 190 individuals was 
generated for the QTL mapping. In brief, a female marine stickleback 
(F 0 ) collected from the Baltic Sea (Helsinki; 60°12'N, 25°11'E) was 
crossed with a male freshwater stickleback (F 0 ) from Lake Pulmanki 
(Lapland; 69°58'N, 27°58'E) in Finland. The female originated from 
a population that most likely represents the ancestral northern Euro- 
pean marine form (Leinonen et al. 201 la). The detailed procedures for 
crossing F 0 grandparents have been described in Leder et al. (2010). 



The resulting F x offspring were reared to maturity in the aquaculture 
facilities of the University of Helsinki. The F x progeny were first raised 
at 17° for 3 mo and then transferred to 4° for 5 mo to simulate 
overwintering, and then transferred back to 17° to stimulate breeding. 
One ¥ l female and one ¥ l male were further crossed to obtain F 2 
progeny. The same F x couple was crossed naturally five times, thus the 
five broods (brood 1: 69 fish; brood 2: 22 fish; brood 3: 22 fish; brood 
4: 44 fish; brood 5: 33 fish) amounted to a total of 190 offspring. Each 
brood was divided into two to six blocks, with an average of 1 1 fish 
per block. Each of the 17 blocks was held in a 27-L aquarium. The F 2 
progeny were reared at 17° for 3 mo, fed ad libitum with frozen 
chironomid larvae, and then killed with an overdose of MS-222. All 
the samples used in the following analyses were preserved in 96% 
ethanol and stored horizontally to avoid bending. The samples were 
then fixed and stained with alizarin red solution for phenotyping. Sex 
of the F 2 offspring was identified by two sex- specific molecular 
markers (GAest31 and Stnl90) based on Natri et al. (2013). 

Morphometric data collection 

Fish preparation, image acquisition, and morphometric measurements 
followed the procedures described in Leinonen et al. (2006). Tradi- 
tional and landmark-based geometric morphometries methods were 
used to quantify variation in body size, body shape, and armor (Figure 
1). In brief, a total of 17 landmarks (for the details, see Leinonen et al. 
2011a) were digitized on the left side of each sample using tpsDig 
version 2.10 (Rohlf 2006), and the aligned coordinate values (X and Y) 
were recorded. Because X and Y landmark coordinates usually map to 
different locations on linkage groups (Albert et al. 2008; Rogers et al. 
2012), they were treated as separate variables. Values of the centroid 
size (Csize) also were recorded as a measure of body size (Bookstein 
1986), and it was strongly correlated with standard length (r s = 0.99, 
n = 185, P < 0.001). Csize is the square root of the sum of the squared 
distances from the measured landmarks to their centroid. This mea- 
sure is independent of any potential random measurement error in the 
landmarks, being a very robust measure of geometric size (Mitteroecker 
and Gunz 2009). In addition, the number of lateral plates on both sides 
(Nplate; the total number of lateral plates on myomeres 1 — 33) was 
quantified from photographs, and four metric variables were measured 
with a digital caliper: (1) length of the first dorsal spine (Dlst), (2) 
length of the second dorsal spine (D2nd), (3) average length of the left 
and right pelvic spines (Pspi), and (4) pelvic girdle length (Pgir). There- 
fore, a total of 40 morphometric variables (six metric and meristic traits, 
including Csize and 34 landmark coordinates) were used in the follow- 
ing QTL mapping analyses. 

DNA extraction and microsatellite genotyping 

Genomic DNA was extracted from fin clips using a silica-based 
purification method (Elphinstone et al. 2003) after proteinase K di- 
gestion. A total of 131 microsatellite markers (Table S2), previously 
isolated for three-spined sticklebacks (Largiader et al. 1999; Peichel 
et al. 2001; Heckel et al. 2002; Colosimo et al. 2004; Colosimo et al. 
2005; Miller et al. 2007; Makinen et al. 2008), were genotyped for the 
grandfather, two ¥ l parents and the 190 F 2 progeny. Genotype data of 
the grandmother (F 0 ) was inferred from the genotypes of twenty F x 
offspring due to sample degradation. Note that primers of all the 
"GAest" markers except GAest46 were from Makinen et al. (2008). 
Primers for GAest46 were forward (5'-AGT GAT CAA TAA CCA 
GAA GGA G-3') and reverse (5'-CGA TAT GCT TTC ATT GTA 
TTT G-3'; H. Makinen, unpublished data). Polymerase chain reac- 
tions (PCRs) were conducted in a 10- uX volume consisting of lx 
QIAGEN Multiplex PCR Master Mix (QIAGEN), 0.5x Q-Solution, 
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Figure 1 Landmark locations for body 
shape and metric and meristic traits. 
Dotted outlines depict the lateral 
plates present in marine populations 
but reduced in freshwater populations. 
The traits and landmarks are described 
in the main text. Nplate, number of 
lateral plates; Pgir, pelvic girdle length; 
Pspi, pelvic spine length; D1st, length 
of the first dorsal spine; D2nd, length of 
the second dorsal spine. 



2 pmol of each primer, and ca. 20 ng of template DNA. The forward 
primers were labeled with FAM, HEX, or TET fluorescent dye, and 
a GTTT-tail was added to the 5 '-end of the reverse primers to pro- 
mote adenylation (Brownstein 1996). PCR conditions were as follows: 
initial denaturation at 95° for 15 min, followed by 30 cycles of 30 sec 
at 94°, 90 sec at 53° and 60 sec at 72°, and a final extension at 60° for 
5 min. PCR products were resolved using a MegaBace 1000 capil- 
lary sequencer (Amersham Biosciences) with ET-ROX 550 size 
standard (Amersham Biosciences) and were analyzed using Frag- 
ment Profiler 1.2 (Amersham Biosciences). All the makers were 
checked for null alleles by using Micro-Checker v.2.2.3 (Van Oosterhout 
et al. 2004). 

Linkage analysis and genetic linkage map construction 

A genetic linkage map was constructed using JoinMap v.4.0 software 
(Van Ooijen 2006). The cross-pollinator population type, which 
allows for segregation of up to four alleles per locus, was used. 
Linkage phases of the loci were determined automatically by the 
software. Test for locus segregation distortion was implemented by 
the x 2 test. All the microsatellites were assigned to linkage groups 
with a relatively stringent two-point logarithm of odds (LOD) score > 
4.0. A map for each linkage group was created using a regression 
mapping module with the following parameters: a recombination 
frequency < 0.499, and a LOD score > 3.0. Map distances (centi- 
Morgans, cM) were calculated using the Kosambi mapping function. 
All the other parameters were set to default values. 

BLAST searches 

To map the physical locations of the microsatellites in the genome, 
BLAST searches were performed against a three-spined stickleback 
genome assembly of Roesti et al. (2013) by using the BLAST module 
within BioEdit v. 7. 1.8 (Hall 1999). This genome assembly improved 
version of the Broad SI assembly of G. aculeatus genome in Ensembl 
Genome Browser (Roesti et al. 2013; http://datadryad.org/resource/ 
doi:10.5061/dryad.846nj). BLAST hits were obtained using the BLASTN 
searching tool with default settings. The expectation value lxe 100 was 
first used to get a unique hit. When no hit was found for a sequence, 
the expectation value decreased to lxe 50 . All the graphic maps were 
drawn using MapChart v.2.2 software (Voorrips 2002). 

Phenotypic data and sexual dimorphism 

A summary of the descriptive statistics for the metric traits is provided 
in Table S3. Because sexual dimorphism was bound to affect the 
morphological variation in our family {e.g., Albert et al. 2008; Leinonen 
et al. 2011b), we tested the effect of sex on body shape (for the land- 
mark coordinates) by means of a multivariate analysis of variance. 
Here, sex, block, and brood were fitted as fixed effects, and Csize as 



covariate. Separate generalized linear mixed models (GLMMs) for 
each metric trait also were run with sex and block as fixed factors, 
Csize as covariate and brood as random effect. GLMM for plate 
counts was run by using a Poisson instead of normal distribution. 
For Csize, GLMM was performed with sex and block as fixed 
factors and brood as a random effect. The analyses were run with 
the multivariate analysis of variance and LMER procedures using 
the statistical software R (R Development Core Team 2008). The 
results revealed that sex had significant effect (P < 0.05) on Csize 
and all the landmark coordinates. 

QTL mapping 

A genome-wide scan for QTL was implemented within R/qtl software 
(Broman et al. 2003). The genotypes, phenotypes, and genetic maps (File 
SI, File S2, and File S3, respectively) converted to four- way cross format 
were imported with the readxross function. Missing genotypes were first 
imputed with fillgeno function. We used the Haley- Knott regression 
method (Haley and Knott 1992) to identify QTL. One-dimensional 
QTL scan (scanone function) was initially carried out to detect QTL 
for landmark coordinates with sex and Csize as covariates. Csize was 
mapped with sex as a covariate. Although sex did not show significant 
effect for the metric and meristic traits, it was also included as a covariate 
for these five traits. Additional QTL were then scanned for each trait 
using a two-QTL model (scantwo function). Finally, a multiple QTL 
model was fitted to determine QTL for each trait using makeqtl and 
fitqtl functions. Locations of QTL were refined with refineqtl function. 
Percentage of phenotypic variance explained (PVE) by each QTL was 
calculated with fitqtl function. The 5% genome- wide significance thresh- 
old was estimated by 1000 permutations. QTL were considered to be 
significant when the LOD scores exceed the 5% genome-wide threshold. 
Confidence interval (CI) for each QTL peak was derived from the 
Bayesian 95% credible interval using the bayesint function. 

RESULTS 

Segregation analysis 

The F 2 progeny comprised 190 offspring, 111 of which were females 
and 79 males. Null alleles were detected in two loci, Stn51 and 
GAest35. Of the 131 loci, 23 (17.6%) showed significantly distorted 
segregation ratios (P < 0.05; Figure 2). Most of them were scattered 
across the linkage groups. However, two clusters of distorted loci were 
found on LG9 (six loci) and LG21 (four loci). Interestingly, most of 
the distorted loci showed missing alleles in the freshwater grandfather. 
Despite the potential effect of distorted segregation ratios on genetic 
linkage map construction and subsequent QTL mapping, all the 
markers were used in the following analyses because the mapping 
analysis was implemented in JoinMap using the independence LOD, 
which is not affected by segregation distortion (Van Ooijen 2006). 
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Genetic linkage map 

The sex- averaged genetic linkage map is shown in Figure 2. Of the 131 
microsatellites, 120 loci were assembled into 21 linkage groups, being 
consistent with earlier karyotypic (Chen and Reisman 1970; Urton 
et al. 2011) and genetic (e.g., Albert et al. 2008; Rogers et al. 2012) 
studies. Of the 1 1 unmapped loci, six appeared to be unlinked to any 
other marker with two-point LOD scores > 4.0 and the other five that 
were assigned to one linkage group were not mapped to any linkage 
group with LOD scores > 3.0 (see details in Table S2). The orientation 
of each linkage group was determined by anchoring at least two loci in 
the published linkage maps of G. aculeatus (e.g., Peichel et al. 2001; 
Albert et al. 2008). The sizes of the linkage groups ranged from 6.4 cM 
(LG15) to 64.4 cM (LG1), spanning in total 652.5 cM of the three- 
spined stickleback genome. The number of markers on each linkage 
group varied between two (LG15) and 12 (LG4) and the average 
intermarker interval was 6.6 cM. 

Comparative genomic mapping 

A total of 129 loci were assembled onto the 21 chromosomes of 
G. aculeatus with extremely low expectation (E)-values < lxe~ 50 
(Figure SI and Table S2). The physical assignment of loci to chro- 
mosomes was consistent with that by genetic linkage analysis with 
four exceptions (viz. Stn34, GAest47, GAest63, and GAest71; Table 
S2). Comparison of the order of loci on the physical map with that 
on the genetic linkage map revealed discordance with at least one 
locus disordered in eight linkage groups/chromosomes (Figure SI). 

QTL mapping 

Results of the genetic linkage groups with significant QTL, 95% CI 
and PVE by QTL for the traits are shown in Figure 2 and Table 1. In 
total, 52 significant QTL were identified. The number of QTL for each 
trait ranged from one to four. A dense set of QTL (23) were mapped 
on LG21, and most of the 95% CIs of QTL overlapped. 

Two significant QTL were detected for Csize, one QTL on LG21 
explaining 27.9% of the variance, and the other on LG6 with PVE 
value of 16.5%. A total of 13 QTL were detected for the armor traits 
(Nplate 3, Dlst 4, D2nd 2, Pgir 2, and Pspi 2; Table 1). One major 
QTL with a rather high PVE (34.6-74.4%) and one or several addi- 
tional QTL with small PVEs (6.0-11.6%) were identified for each of 
these traits (Table 1). For Nplate, one QTL mapped to locus close to 
Stn380 on LG4, showed the highest PVE value (74.4%), and the other 
two QTL had minor effect (9.1-10.6%). Major QTL for the other four 
armor traits all mapped on LG21 with overlapping 95% CIs, and the 
PVEs ranged from 34.6 (Pgir) to 57.6% (D2nd). 

The X and Y coordinates of each landmark usually mapped to 
different linkage group regions, implying their distinctiveness in the 
analyses. A total of 37 significant QTL were found for the 34 X/Y 
coordinates, 17 of which were mapped to LG21. Of the 37 QTL, 21 
were detected for X-coordinates and 16 for Y-coordinates. Ten QTL 
were found with major effect (PVE > 15%), which were all located on 
LG21. Eight QTL with major effect were found for X-coordinates. 

DISCUSSION 

In the present study, we discovered multiple QTL in the three-spined 
stickleback, some of which had large phenotypic effects, and many 



more had moderate- to- small effects. The findings appear to be, to 
some extent, consistent with those obtained from earlier Pacific 
crosses of this species, but also many novel QTL were discovered. In 
what follows, we discuss these issues in more detail and relate our 
findings to those that have emerged from earlier QTL studies of three- 
spined sticklebacks. 

QTL mapping of the ecologically important traits 

Body size, body shape, and armor structures are ecologically im- 
portant traits frequently under directional natural and/or sexual 
selection (Walker 1997; Blanckenhorn 2000; Andersson et al. 2006; 
Head et al. 2009). These phenotypic traits show genetically based 
differentiation among different three-spined stickleback populations 
(e.g., Wright et al. 2004; McGuigan et al. 2011; Karhunen et al. 2013). 
During the last decade, studies have started to uncover the genetic bases 
underlying these traits, especially the armor traits such as the lateral 
plates and the pelvic structure (Peichel et al. 2001; Colosimo et al. 2004; 
Cresko et al. 2004; Shapiro et al. 2004; Coyle et al. 2007). In particular, 
two major genes — Ectodysplasin (Eda) and Pituitary homebox transcrip- 
tion factor 1 (Pitxl) — have been identified to account for plate and 
pelvic reduction in three-spined sticklebacks, respectively (Shapiro 
et al. 2004; Colosimo et al. 2005; Coyle et al. 2007; Chan et al. 2010). 

In this study, 52 significant QTL were identified for traits con- 
cerning the body size (centroid size), body shape (X and Y landmark 
coordinates), and body armor. The pattern of one major plus 
several minor QTL was observed for most traits, as is the case also 
in previous studies of sticklebacks (e.g., Colosimo et al. 2004; Shapiro 
et al. 2004). One major QTL on LG4 was detected for lateral plate 
number, with the LOD peak close to Stn380, a marker located within 
Eda gene. This finding is consistent with one earlier study, where one 
major QTL was also found on LG4 with LOD peak close to marker 
Gac4174 (Colosimo et al. 2004) near the Eda gene (see Table S2 for 
physical locations of Stn380 and Gac4174). However, major QTL for 
body size, body shape and armor traits (excluding Nplate) were all 
mapped to LG21, being different from earlier mapping results (Table 
2). Such differences also have been observed in previous studies (Table 
2). In addition, modifier QTL with smaller effect detected in this and 
some previous studies also seem to vary between different stickleback 
populations (Table 2). These differences may be caused by their geo- 
graphically different parental origins and possibly also by the narrow 
representation of allelic variation in crosses based on a small number 
of individuals. Nevertheless, the results also are consistent with the 
possibility that different genetic architectures may underlie expression 
of morphological traits in stickleback populations of different geo- 
graphic origins. 

As for the QTL effect sizes, their distributions appear to be quite 
comparable in different stickleback crosses, although the QTL are 
mapped often to different linkage groups (Table 2; Figure S2). For 
body size (i.e., centroid size), the major QTL explained -30% of the 
variance in this study, and for its other proxy, the standard body 
length, the major QTL explained 21% and 41% of the variance in 
two different Pacific crosses (Kitano et al. 2009; Greenwood et al. 
2011). The major QTL for lateral plate number explained more than 
70% of the variance in this and one previous study (Colosimo et al. 
2004). For body shape, several major QTL had PVEs up to 40% in this 



Figure 2 A genome-wide consensus linkage map of G. aculeatus and the significant quantitative trait loci detected for phenotypic traits. Traits are 
divided into body size (centroid size), body shape (landmark coordinates, X/Y), and armor traits with different colors. See Figure 1 for trait 
abbreviations. The locations of two genes, Eda for lateral plate number and Pitxl for pelvic development, are shown in orange. Eda includes 
Stn380, and the location of Pitxl is below Stn82. Asterisks indicate markers with significantly distorted segregation ratios (*P < 0.05; **P < 0.01). 
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Table 1 Significant QTL for the meristic, metric traits, and landmark coordinates in three-spine sticklebacks 



Trait 


LG 


Location, cM 


Nearest locus, cM 


LOD 


PVE, % 


95% CI, cM 


Nplate 


4 


12.9 


Stn380 (12.9) 


71.9 


74.4 


12.9-13.0 




9 


12.0 


Stn108 (9.8) 


18.9 


9.1 


12.0-12.0 




21 


8.0 


Stn223 (7.2) 


21.3 


10.6 


7.0-10.0 


D1st 


7 


26 


Stn76 (26.8) 


8.72 


9.05 


23.0-26.0 




11 


13.2 


Gac7033 (13.2) 


10.9 


11.64 


13.0-25.0 




16 


0 


Stn315 (0.0) 


6.0 


6.01 


0.0-15.0 




21 


8 


Stn223 (7.2) 


34.2 


50.55 


8.0-9.0 


D2nd 


7 


47.0 


Stn73 (47.2) 


6.75 


7.26 


1.0-53.0 




21 


8.0 


Stn223 (7.2) 


36.32 


57.64 


8.0-10.0 


Pspi 


7 


49.0 


Stn73 (47.2) 


5.86 


6.32 


1.0-55.0 




21 


9.0 


Stn223 (7.2) 


35.29 


56.16 


8.0-10.0 


Pgir 


13 


11.3 


Stn154 (11.3) 


6.96 


10.44 


0.0-9.0 




21 


11.0 


GAest8 (1 1 .0) 


19.62 


34.60 


8.0-11.0 


Csize 


6 


9.0 


Stn306 (8.7) 


8.92 


16.53 


5.0-28.0 




21 


9.0 


Stn223 (7.2) 


14.08 


27.89 


8.0-10.0 


X1 


13 


11.3 


Stn154 (11.3) 


4.69 


9.32 


4.0-11.3 




21 


9.0 


Stn223 (7.2) 


8.10 


16.81 


6.0-11.0 


Y1 


21 


11.0 


GAest8 (1 1 .0) 


4.84 


7.60 


6.0-11.0 


X2 


1 


10.6 


Stn5 (10.6) 


4.88 


8.53 


5.0-21.0 




10 


15.2 


GAest46 (15.2) 


4.88 


8.53 


14.0-25.0 




21 


10.0 


GAest8 (1 1 .0) 


10.19 


19.06 


8.0-11.0 


Y3 


13 


11.3 


Stn154 (11.3) 


4.61 


8.0 


6.0-18.0 




17 


36.0 


GAest4 (37.3) 


4.46 


7.72 


26.0-37.3 




21 


9.0 


Stn223 (7.2) 


4.48 


7.78 


6.0-11.0 


X4 


13 


17.0 


Stn156 (16.5) 


6.50 


8.92 


0.0-18.0 




15 


1.0 


Stn170 (0.0) 


4.36 


5.83 


0.0-5.0 




20 


8.0 


Stn214 (5.9) 


6.00 


8.18 


5.9-9.0 




21 


8.0 


Stn223 (7.2) 


16.50 


25.82 


8.0-10.0 


Y4 


16 


4.0 


Stn174 (3.4) 


4.58 


8.06 


0.0-16.0 




21 


8.0 


Stn223 (7.2) 


8.46 


15.64 


6.0-11.0 


X5 


13 


11.0 


Stn154 (11.3) 


5.42 


7.09 


2.0-17.0 




21 


8.0 


Stn223 (7.2) 


24.41 


41.05 


8.0-9.0 


Y5 


12 


27.4 


GAest8 (1 1 .0) 


5.82 


11.66 


25.8-54.0 




16 


5.0 


Stn174 (3.4) 


5.07 


10.05 


3.0-7.3 


X6 


21 


8.0 


Stn223 (7.2) 


6.26 


12.81 


8.0-10.0 


Y6 


5 


20.0 


Stn52 (18.1) 


5.25 


8.85 


19.0-25.0 




12 


31.0 


Stn146 (32.7) 


5.38 


9.10 


27.4-34.0 




19 


12.0 


Stn194 (12.1) 


4.72 


7.90 


0.0-13.0 


X8 


21 


11.0 


GAest8 (1 1 .0) 


11.95 


24.16 


6.0-11.0 


Y8 


13 


0.0 


GAest67 (0.0) 


4.42 


9.89 


1.0-35.0 


X9 


21 


9.0 


Stn223 (7.2) 


8.17 


15.07 


8.0-11.0 


X11 


10 


20.0 


GAest63 (15.8) 


5.68 


11.52 


16.0-33.0 




12 


36.0 


Stn327 (38.2) 


4.67 


9.34 


35.0-37.0 


Y11 


21 


9.0 


Stn223 (7.2) 


16.69 


32.53 


8.0-10.0 


X12 


1 


26.0 


Stn8 (26.8) 


4.54 


7.74 


19.0-33.0 




21 


9.0 


Stn223 (7.2) 


12.23 


23.51 


7.0-11.0 


Y12 


13 


0.0 


GAest67 (0.0) 


4.79 


10.02 


0.0-5.0 




21 


8.0 


Stn223 (7.2) 


5.97 


12.67 


7.2-10.0 


X13 


21 


9.0 


Stn223 (7.2) 


18.41 


34.47 


8.0-10.0 


Y15 


21 


11.0 


GAest8 (1 1 .0) 


6.24 


13.07 


8.0-11.0 


X16 


21 


11.0 


GAest8 (1 1 .0) 


5.96 


12.70 


7.0-11.0 


X17 


21 


3.0 


Stn208 (3.5) 


5.50 


11.89 


0.0-9.0 



QTL, quantitative trait locus; LG, linkage group; LOD, logarithm of odds; PVE, percentage of phenotypic variance explained; 95% CI, 95% confidence interval. See 
Figure 1 for trait abbreviations 



and two previous studies (Albert et al. 2008; Rogers et al. 2012). In 
addition, it seems that the minor QTL for these morphological traits 
mostly explained -5-10% of the variance. QTL with PVE below 5% 
are rare, probably because QTL with small effect size are difficult to 
detect with genetic mapping approaches such as the one utilized here. 
Furthermore, it is interesting to observe that the distributions of QTL 
effect sizes (PVEs) for shape traits between different crosses/studies 
are similar (Figure S2). This observation supports the conjecture that 



different chromosomal regions might influence the body shape in 
different populations, but the general genetic architecture in terms 
of QTL effect sizes is very similar. 

Multiple QTL on LG21 

Among the 21 linkage groups, quite a few QTL (44%) for body size, 
armor and body shape were mapped on LG21, and many of them 
appeared to colocalize {i.e., the peak LOD in the same location and the 
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Table 2 Comparison of significant QTL for morphological traits used here between this and the previous studies of three-spined 
sticklebacks 



Trait 




This Study 








Previous Study 




KJ. I L 


i 


p\/p 0/ 

rvt, /o 


rceT 


\J. I L 


I 


p\/p 0/ 

rvt, /o 


Body size 
















Csize/standard 


Major 


LG21 


27.9 


2 


Major 


LG19 


13 


body length 


Minor 


LG6 


10.4 


5 




LG13 


NA 










6 


Major 


LG19 


20.7 












Minor 


LG19 


12.1 










7 


Major 


LG17, LG18 


23.6, 41.1 


Armor trait 
















Nplate 


Major 


LG4 


74.4 


1 


Major 


LG13, LG21 b 


26, 23 




Minor 


LG9, LG21 


9.1-10.6 


2 


Major 


LG4 


76.9 












Minor 


LG7, LG10, LG21 


3.7-6.5 


D1st 


Major 


LG21 


50.6 


1 


Major 


LG1 , LG2 


21, 17 




Minor 


LG7, LG11, LG16 


6.0-10.9 


6 


Major 


LG9 


45.8, 50.4 


D2nd 


Major 


LG21 


57.6 


1 


Major 


LG8, LG11 


22, 17 




Minor 


LG7 


7.3 










Pspi 


Major 


LG21 


56.2 


1 


Major 


LG8 


25 




Minor 


LG7 


6.3 


3 


Major 


LG7 


65.3, 43.7 












Minor 


LG2, LG4 


7.6, 5.8 










4 


Major 


LG7 


85.5 


Pgir 


Major 


LG21 


34.6 


3 


Major 


LG7 


46.8, 27.8 




Minor 


LG13 


10.4 




Minor 


LG1 , LG2, LG4 


5.6-11.1 










4 


Major 


LG7 


87 


Body shape 
















Coordinate 


Major 


LG21 


15.1-41.1 


5 


Major 


LG19, LG7, LG4 


15.8-36 


traits 


















Minor 


LG1, LG5, LG10, LG12, 


5.8-13.1 




Minor 


LG1-5, LG7-LG9, LG12, 


2.5-14.6 






LG13, LG5-LG17, 








LG13, LG15-LG21 








LG19-LG21 




8 


Major 


LG4 


15.1-40.6 












Minor 


LG1 , LG2, LG4, LG7-LG14, 


5.8-14.7 



LG16-LG21 



Major QTL are defined as PVE > 15%, whereas minor QTL are PVE < 15%. QTL, quantitative trait locus; PVE, percentage of phenotypic variance explained; NA, not 
available. 

a Ref, references including 1: Peichel etal. (2001); 2: Colosimo etal. (2004); 3: Shapiro etal. (2004); 4: Coyle etal. (2007); 5: Albert etal. (2008); 6: Kitano etal. (2009); 
b 7: Greenwood et a/. (201 1); and 8: Rogers et al. (2012). 
For comparison, the same linkage group names were used. 



95% CIs overlapping). It is probably not a coincidence that most of 
these traits were also significantly intercorrelated at phenotypic level 
(Table S4). Either pleiotropic effects of a single QTL or tight linkage 
between genes that influence different morphological traits might 
contribute to the clustering. Irrespective of the actual cause, the result 
suggests that the QTL on LG21 may have an important role in con- 
trolling many adaptive traits in three-spined sticklebacks. Interestingly, 
two markers (Stn208, 6.5 Mb; Stn223, 7.5 Mb) on LG21 located close to 
peaks of these overlapping QTL were in a region where an inversion 
(5.8-7.5 Mb) on chromosome XXI (i.e., LG21) between Paxton benthic 
and Japanese Pacific sticklebacks has been recently reported (Jones et al. 
2012). Several QTL controlling lateral plate numbers (Colosimo et al. 
2004), body shape traits (Albert et al. 2008), and lateral line-related 
traits (Wark et al. 2012) also have been found to map to this same 
inversion. Given the low marker density on LG21, more markers should 
be developed and fine-scale mapping should be conducted to better 
understand the genetic structure of this QTL hotspot. 

QTL estimation bias 

Potential biases in the interpretation of QTL effect sizes may be 
attributable to, for instance, the Beavis effect, referring to over- 
estimation of QTL effects due to small sample sizes (n < 500; Beavis 
1994, 1998). The Beavis effect is caused by the fact that QTL are 



reported when the test statistics reach a predetermined threshold, 
and thus, their estimated effects tend to be upward-biased (Barton 
and Keightley 2002; Xu 2003). Therefore, increasing the number of 
individuals and/or markers is expected to detect QTL with smaller 
effect size (Mackay et al. 2009). The Beavis effect can be alleviated by 
replicating experiments, or by estimating QTL effect sizes from a dif- 
ferent sample of individuals instead of that used for QTL detection 
(Slate 2005). 

In three-spined sticklebacks, QTL for similar or identical morpho- 
logical traits have been identified in some earlier crosses. For instance, 
comparable distributions of QTL effect size for body shape traits in this 
and earlier studies (Albert et al. 2008; Rogers et al. 2012; Table 2 and 
Figure S2) suggest that the influence of Beavis effect on body shape 
traits in this study has probably been small. Furthermore, QTL effect 
sizes for lateral plate number estimated in this study were very similar 
to those detected in an earlier study (Colosimo et al. 2004; Table 2). 
Therefore, it appears that the Beavis effect might have little influence on 
the QTL detection power in this study. 

Linkage map vs. physical map 

Four markers were assigned to linkage groups that differed from the 
chromosomes in this study. And several conflicts also were detected in 
marker order between the linkage and physical maps. Several reasons 
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might explain these discrepancies. The first is the accuracy of the 
linkage mapping approach, which can be affected by factors such as 
genotyping errors, missing data and/or segregation distortion (Hackett 
and Broadfoot 2003; Slate 2008). For example, Stn34 was mapped on 
LG3 in Peichel et al. (2001), and on chromosome III in our BLAST 
searches, but on LG4 in our genetic linkage map. Second, there may be 
errors in the reference genome sequence (Ross and Peichel 2008; Natri 
et al. 2013; Roesti et al. 2013). Third, the disagreements might stem 
from complex recombination events between divergent genomes, 
rather than from technical errors. However, some of the loci showing 
discrepancies {e.g., GAest47, GAest63, and GAest71) were used in the 
genetic linkage mapping for the first time, and thus, comparisons to 
other crosses are not possible. 

Distorted segregation ratios of microsatellite markers 

The phenomenon that the segregation ratio of a locus deviates from 
the expected Mendelian ratio is of common occurrence in genome 
analyses (Lu et al 2002; Xu 2008). Technical problems such as null 
alleles can yield distorted marker segregation ratios. Here, null alleles 
were detected in two microsatellite loci (Stn51 and GAest35). As ex- 
pected, these two loci showed extremely distorted segregation ratios 
(P < 0.01). Alternatively, a variety of genetic and/or physiological 
factors, including differential transmission in male or female germ 
line and post-zygotic selection prior to genotypic evaluation (Xu 
et al 1997), can yield distorted segregation ratios as well. Genetic 
incompatibility is one of the most common reasons, and increased 
divergence between species or populations is expected to increase 
the deviations from the expected Mendelian segregation (Zamir and 
Tadmor 1986). Although the marine and freshwater G. aculeatus 
populations are phenotypically and genetically divergent, their off- 
spring are highly viable and fertile, as shown in this and earlier studies 
(e.g., Peichel et al. 2001; Colosimo et al. 2004; Shapiro et al. 2004; 
Albert et al. 2008; Rogers et al. 2012). Therefore, genetic incompati- 
bilities seem to be an unlikely explanation in our case. 

Presence of segregation distortion loci (SDL; Vogl and Xu 2000) 
also could lead to distorted segregation ratios of molecular markers. 
SDL are loci that are subject to gametic or zygotic selection and 
consequently cause segregation distortion (Xu 2008). Loci linked to 
SDL can be indirectly distorted because of genetic hitchhiking (Xu 
et al. 1997), thus leading to clustering of distorted loci. This phenom- 
enon is detected on LG9 and LG21 in our linkage map. The clustering 
of the distorted markers suggests that deleterious recessive alleles are 
probably linked to them (Zamir and Tadmor 1986). Regions with 
three or more closely linked loci that exhibit significant segregation 
distortion are termed segregation distortion regions (SDRs), and iden- 
tifying common SDRs among different populations would be helpful 
to determine the underlying SDL (Lu et al. 2002). However, we are not 
certain that the two regions are actually common SDRs among dif- 
ferent three-spined stickleback crosses, because no distorted markers 
have ever been reported in the genetic linkage maps of the three- 
spined stickleback (e.g., Peichel et al. 2001; Albert et al. 2008; Rogers 
et al. 2012). 

In summary, we identified multiple QTL for a set of ecologically 
important morphological traits that are known to exhibit significant 
differentiation between marine and freshwater stickleback popula- 
tions. The results identified several significant QTL with large effects 
in several morphological traits, as well as a large number of QTL with 
smaller effects. In accordance with earlier studies, a QTL on LG4 was 
identified to have a major effect on the lateral plate number. In 
addition, several major overlapping QTL were mapped on LG21, 
suggesting pleiotropic effects or tight genetic linkage between QTL. 



These QTL were associated with centroid size, lengths of pelvic girdle, 
pelvic and dorsal spines, and different aspects of body shape, implying 
the potential role of genomic segments on LG21 in multiple adaptive 
processes. Further studies are needed to identify candidate genes on 
this linkage group (i.e., chromosome 21) controlling the correlated 
complex morphological traits of three-spined stickleback. All in all, 
this study provides new insights into genetic architecture of adaptive 
morphological traits in three-spined sticklebacks and should aid fur- 
ther studies aimed at deciphering the genetic basis of morphological 
variation in this species. 

ACKNOWLEDGMENTS 

We thank Marika Karjalainen and Linda Uoti for their help at various 
stages of this work and Jacquelin DeFaveri, the editor, and two 
anonymous reviewers for comments that improved the manuscript. 
This work was supported by the Chinese Academy of Sciences (the 
100-talent Program of the Chinese Academy of Sciences; M.-H.L.), 
the Academy of Finland (No. 250435 and No. 265211; J.M.) and the 
China Postdoctoral Science Foundation (No. 055S004B04; J.L.). 

LITERATURE CITED 

Albert, A. Y., S. Sawaya, T. H. Vines, A. K. Knecht, C. T. Miller et al, 
2008 The genetics of adaptive shape shift in stickleback: pleiotropy and 
effect size. Evolution 62: 76-85. 

Andersson, J., F. Johansson, and T. Soderlund, 2006 Interactions between 
predator- and diet-induced phenotypic changes in body shape of crucian 
carp. Proc. R. Soc. B 273: 431-437. 

Barrett, R. D., 2010 Adaptive evolution of lateral plates in three-spined 
stickleback Gasterosteus aculeatus: a case study in functional analysis of 
natural variation. J. Fish Biol. 77: 311-328. 

Barton, N. H., and P. D. Keightley, 2002 Understanding quantitative 
genetic variation. Nat. Rev. Genet. 3: 11-21. 

Beavis, W. D., 1994 The power and deceit of QTL experiments: lessons 
from comparative QTL studies, pp. 250-266 in Proceedings of the Forty- 
Ninth Annual Corn and Sorghum Industry Research Conference, edited by 
D. B. Wilkinson. American Seed Trade Association, Chicago, IL. 

Beavis, W. D., 1998 QTL analyses: power, precision and accuracy, pp. 145- 
162 in Molecular Dissection of Complex Traits, edited by A. H. Paterson, 
C. R. C. Press, Boca Raton, FL. 

Bell, M. A., and S. A. Foster, 1994 The Evolutionary Biology of the 
Threespine Stickleback, Oxford University Press, Oxford. 

Bergelson, J., and F. Roux, 2010 Towards identifying genes underlying 
ecologically relevant traits in Arabidopsis thaliana. Nat. Rev. Genet. 11: 
867-879. 

Blanckenhorn, W., 2000 The evolution of body size: what keeps organisms 

small? Q. Rev. Biol. 75: 385-407. 
Bookstein, F., 1986 Size and shape spaces for landmark data in two 

dimensions. Stat. Sci. 1: 181-242. 
Broman, K. W., H. Wu, S. Sen, and G. A. Churchill, 2003 R/qtl: QTL 

mapping in experimental crosses. Bioinformatics 19: 889-890. 
Brownstein, M. J., 1996 Modulation of non-templated nucleotide addition 

by Taq DNA polymerase: primer modifications that facilitate genotyping. 

Biotechniques 20: 1004-1010. 
Chan, Y. F., M. E. Marks, F. C. Jones, G. Villarreal, Jr, M. D. Shapiro et al, 

2010 Adaptive evolution of pelvic reduction in sticklebacks by recurrent 

deletion of a Pitxl enhancer. Science 327: 302-305. 
Chen, C, and K. Ritland, 2013 Lineage-specific mapping of quantitative 

trait loci. Heredity 111: 106-113. 
Chen, T. R., and H. M. Reisman, 1970 A comparative chromosome study of 

the North American species of sticklebacks (Teleostei: Gasterosteidae). 

Cytogenet. Genome Res. 9: 321-332. 
Colosimo, P. F., K. E. Hosemann, S. Balabhadra, G. Villarreal, Jr, M. Dickson 

et al, 2005 Widespread parallel evolution in sticklebacks by repeated 

fixation of Ectodysplasin alleles. Science 307: 1928-1933. 



602 I J. Liu eta/. 



:~-ZLG3' Genes | Genomes | Genetics 



Colosimo, P. F., C. L. Peichel, K. Nereng, B. K. Blackman, M. D. Shapiro 

et al, 2004 The genetic architecture of parallel armor plate reduction in 

threespine sticklebacks. PLoS Biol. 2: El 09. 
Coyle, S. M., F. A. Huntingford, and C. L. Peichel, 2007 Parallel evolution of 

Pitxl underlies pelvic reduction in Scottish threespine stickleback 

(Gasterosteus aculeatus). J. Hered. 98: 581-586. 
Cresko, W. A., A. Amores, C. Wilson, J. Murphy, M. Currey et al, 

2004 Parallel genetic basis for repeated evolution of armor loss in 
Alaskan threespine stickleback populations. Proc. Natl. Acad. Sci. USA 
101: 6050-6055. 

Ellegren, H., and B. C. Sheldon, 2008 Genetic basis of fitness differences in 

natural populations. Nature 452: 169-175. 
Elphinstone, M. S., G. N. Hinten, M. J. Anderson, and C. J. Nock, 2003 An 

inexpensive and high-throughput procedure to extract and purify total 

genomic DNA for population studies. Mol. Ecol. Notes 3: 317-320. 
Erickson, D. L., C. B. Fenster, H. K. Stenoien, and D. Price, 2004 Quantitative 

trait locus analyses and the study of evolutionary process. Mol. Ecol. 13: 

2505-2522. 

Gibson, G., 2005 The synthesis and evolution of a supermodel. Science 307: 
1890-1891. 

Greenwood, A. K., F. C. Jones, Y. F. Chan, S. D. Brady, D. M. Absher et al, 

2011 The genetic basis of divergent pigment patterns in juvenile 
threespine sticklebacks. Heredity 107: 155-166. 

Hackett, C. A., and L. B. Broadfoot, 2003 Effects of genotyping errors, 
missing values and segregation distortion in molecular marker data on 
the construction of linkage maps. Heredity 90: 33-38. 

Haley, C. S., and S. A. Knott, 1992 A simple regression method for mapping 
quantitative trait loci in line crosses using flanking markers. Heredity 69: 
315-324. 

Hall, T. A., 1999 BioEdit: a user-friendly biological alignment editor and anal- 
ysis program for Windows 95/98/NT. Nucl. Acids. Symp. Ser. 41: 95-98. 

Head, M. L., E. A. Price, and J. W. Boughman, 2009 Body size differences 
do not arise from divergent mate preferences in a species pair of 
threespine stickleback. Biol. Lett. 5: 517-520. 

Heckel, G., M. Zbinden, D. Mazzi, A. Kohler, G. Reckeweg et al, 

2002 Microsatellite markers for the three-spined stickleback (Gasteros- 
teus aculeatus L.) and their applicability in a freshwater and an anadro- 
mous population. Conserv. Genet. 3: 79-81. 

Jones, F. C., M. G. Grabherr, Y. F. Chan, P. Russell, E. Mauceli et al, 

2012 The genomic basis of adaptive evolution in threespine stickle- 
backs. Nature 484: 55-61. 

Karhunen, M., J. Merila, T. Leinonen, J. M. Cano, and O. Ovaskainen, 

2013 DRIFTSEL: an R package for detecting signals of natural selection 
in quantitative traits. Mol. Ecol. Res. 13: 746-754. 

Karlsson, E. K., and K. Lindblad-Toh, 2008 Leader of the pack: gene 
mapping in dogs and other model organisms. Nat. Rev. Genet. 9: 
713-725. 

Kimmel, C. B., B. Ullmann, C. Walker, C. Wilson, M. Currey et al, 

2005 Evolution and development of facial bone morphology in threes- 
pine sticklebacks. Proc. Natl. Acad. Sci. USA 102: 5791-5796. 

Kitano, J., J. A. Ross, S. Mori, M. Kume, F. C. Jones et al, 2009 A role for 
a neo-sex chromosome in stickleback speciation. Nature 461: 1079-1083. 

Largiader, C. R., V. Fries, B. Kobler, and T. C. Bakker, 1999 Isolation and 
characterization of microsatellite loci from the three-spined stickleback 
(Gasterosteus aculeatus L.). Mol. Ecol. 8: 342-344. 

Lavin, P. A., and J. D. McPhail, 1986 Adaptive divergence of trophic phe- 
notype among freshwater populations of the threespine stickleback 
(Gasterosteus aculeatus). Can. J. Fish. Aquat. Sci. 43: 2456-2463. 

Leder, E. H., J. M. Cano, T. Leinonen, R. B. O'Hara, M. Nikinmaa et al, 
2010 Female-biased expression on the X chromosome as a key step in 
sex chromosome evolution in threespine sticklebacks. Mol. Biol. Evol. 27: 
1495-1503. 

Leinonen, T., J. M. Cano, H. Makinen, and J. Merila, 2006 Contrasting 
patterns of body shape and neutral genetic divergence in marine and lake 
populations of threespine sticklebacks. J. Evol. Biol. 19: 1803-1812. 

Leinonen, T., J. M. Cano, and J. Merila, 2011a Genetics of body shape and 
armour variation in threespine sticklebacks. J. Evol. Biol. 24: 206-218. 



Leinonen, T., J. M. Cano, and J. Merila, 2011b Genetic basis of sexual 
dimorphism in the threespine stickleback Gasterosteus aculeatus. He- 
redity 106: 218-227. 

Lu, H., J. Romero -Severson, and R. Bernardo, 2002 Chromosomal regions 
associated with segregation distortion in maize. Theor. Appl. Genet. 105: 
622-628. 

Lynch, M., and B. Walsh, 1998 Genetics and Analysis of Quantitative Traits. 
Sinauer, Sunderland, MA. 

Mackay, T. F., E. A. Stone, and J. F. Ayroles, 2009 The genetics of quan- 
titative traits: challenges and prospects. Nat. Rev. Genet. 10: 565-577. 

Makinen, H., J. Cano, and J. Merila, 2008 Identifying footprints of directional 
and balancing selection in marine and freshwater three-spined 
stickleback (Gasterosteus aculeatus) populations. Mol. Ecol. 17: 3565-3582. 

McGuigan, K., N. Nishimura, M. Currey, D. Hurwit, and W. A. Cresko, 

2011 Cryptic genetic variation and body size evolution in threespine 
stickleback. Evolution 65: 1203-1211. 

McKinnon, J. S., and H. D. Rundle, 2002 Speciation in nature: the threes- 
pine stickleback model systems. Trends Ecol. Evol. 17: 480-488. 

McKinnon, J. S., S. Mori, B. K. Blackman, L. David, D. M. Kingsley et al, 
2004 Evidence for ecology's role in speciation. Nature 429: 294-298. 

Merila, J., and B. C. Sheldon, 1999 Genetic architecture of fitness and nonfitness 
traits: empirical patterns and development of ideas. Heredity 82: 103-109. 

Miller, C. T., S. Beleza, A. A. Pollen, D. Schluter, R. A. Kittles et al, 

2007 ds-regulatory changes in Kit ligand expression and parallel evo- 
lution of pigmentation in sticklebacks and humans. Cell 131: 1179-1189. 

Mitteroecker, P., and P. Gunz, 2009 Advances in geometric morphometries. 
Evol. Boil. 36: 235-247. 

Mott, R., C. J. Talbot, M. G. Turri, A. C. Collins, and J. Flint, 2000 A 
method for fine mapping quantitative trait loci in outbred animal stocks. 
Proc. Natl. Acad. Sci. USA 97: 12649-12654. 

Natri, H. M., T. Shikano, and J. Merila, 2013 Progressive recombination 
suppression and differentiation in recently evolved neo-sex chromo- 
somes. Mol. Biol. Evol. 30: 1131-1144. 

Nordborg, M., and D. Weigel, 2008 Next-generation genetics in plants. 
Nature 456: 720-723. 

Peichel, C. L., K. S. Nereng, K. A. Ohgi, B. L. E. Cole, P. F. Colosimo et al, 
2001 The genetic architecture of divergence between threespine 
stickleback. Nature 414: 901-905. 

Phillips, P. C, 2005 Testing hypotheses regarding the genetics of 
adaptation. Genetica 123: 15-24. 

Price, A. H., 2006 Believe it or not, QTLs are accurate! Trends Plant Sci. 11: 
213-216. 

Roesti, M., D. Moser, and D. Berner, 2013 Recombination in the threespine 
stickleback genome-patterns and consequences. Mol. Ecol. 22: 3014- 
3027. 

Rogers, S., P. Tamkee, B. Summers, S. Balabahadra, M. Marks et al, 

2012 Genetic signature of adaptive peak shift in threepine stickleback. 
Evolution 66: 2439-2450. 

Rohlf, F. J., 2006 TPS Software Series (distributed by the author). Depart- 
ment of Ecology and Evolution, State University of New York, Stony 
Brook, NY. 

Ross, J. A., and C. L. Peichel, 2008 Molecular cytogenetic evidence of re- 
arrangements on the Y chromosome of the threespine stickleback fish. 
Genetics 179: 2173-2182. 

Shapiro, M. D., M. E. Marks, C. L. Peichel, B. K. Blackman, K. S. Nereng 
et al, 2004 Genetic and developmental basis of evolutionary pelvic re- 
duction in threespine sticklebacks. Nature 428: 717-724. 

Slate, J., 2005 Quantitative trait locus mapping in natural populations: 
progress, caveats and future directions. Mol. Ecol. 14: 363-379. 

Slate, J., 2008 Robustness of linkage maps in natural populations: a simu- 
lation study. Proc. R. Soc. B 275: 695-702. 

Slate, J., 2013 From beavis to beak color: a simulation study to examine how 
much qtl mapping can reveal about the genetic architecture of quanti- 
tative traits. Evolution 67: 1251-1262. 

Slate, J., A. W. Santure, P. G. Feulner, E. A. Brown, A. D. Ball et al, 

2010 Genome mapping in intensively studied wild vertebrate popula- 
tions. Trends Genet. 26: 275-284. 



G3'Genes | Genomes | Genetics 



Volume 4 April 2014 I Three-Spined Stickleback QTL Mapping I 603 



Symonds, V. V., A. V. Godoy, T. Alconada, J. F. Botto, T. E. Juenger et al, 
2005 Mapping quantitative trait loci in multiple populations of Arabi- 
dopsis thaliana identifies natural allelic variation for trichome density. 
Genetics 169: 1649-1658. 

Urton, J. R., S. R. McCann, and C. L. Peichel, 2011 Karyotype differentia- 
tion between two stickleback species (Gasterosteidae). Cytogenet. Ge- 
nome Res. 135: 150-159. 

Van Ooijen, J., 2006 JoinMap 4.0, Software for the Calculation of Genetic 
Linkagemaps in Experimental Populations. Kyazma B. V., Wageningen, 
Netherlands. 

Van Oosterhout, C., W. F. Hutchinson, D. P. M. Wills, and P. Shipley, 
2004 Micro-checker: software for identifying and correcting 
genotyping errors in microsatellite data. Mol. Ecol. Notes 4: 535- 
538. 

Vogl, C., and S. Xu, 2000 Multipoint mapping of viability and segre- 
gation distorting loci using molecular markers. Genetics 155: 1439- 
1447. 

Voorrips, R., 2002 MapChart: software for the graphical presentation of 

linkage maps and QTLs. J. Hered. 93: 77-78. 
Walker, J., 1997 Ecological morphology of lacustrine threespine stickleback 

Gasterosteus aculeatus L. (Gasterosteidae) body shape. Biol. J. Linn. Soc. 

Lond. 61: 3-50. 



Walker, J., and M. A. Bell, 2000 Net evolutionary trajectories of body shape 
evolution within a microgeographic radiation of threespine sticklebacks 
{Gasterosteus aculeatus). J. Zool. (Lond.) 252: 293-302. 

Wark, A., M. Mills, L. H. Dang, Y. F. Chan, F. C. Jones et al, 2012 Genetic 
architecture of variation in the lateral line sensory system of threespine 
sticklebacks. G3 (Bethesda) 2: 1047-1056. 

Wright, H. A., R. J. Wootton, and I. Barber, 2004 Interpopulation variation 
in early growth of threespine sticklebacks (Gasterosteus aculeatus) under 
laboratory conditions. Can. J. Fish. Aquat. Sci. 61: 1832-1838. 

Wu, R., and M. Lin, 2006 Functional mapping — how to map and study the 
genetic architecture of dynamic complex traits. Nat. Rev. Genet. 7: 229-237. 

Xu, S., 2003 Theoretical basis of the Beavis effect. Genetics 165: 2259-2268. 

Xu, S., 2008 Quantitative trait locus mapping can benefit from segregation 
distortion. Genetics 180: 2201-2208. 

Xu, Y., L. Zhu, N. Huang, and S. R. McCouch, 1997 Chromosomal regions 
associated with segregation distortion of molecular markers in F 2 , back- 
cross, doubled haploid and recombinant inbred populations in rice (Or- 
yza satiava L.). Mol. Gen. Genet. 253: 535-545. 

Zamir, D., and Y. Tadmor, 1986 Unequal segregation of nuclear genes in 
plants. Bot. Gaz. 147: 355-358. 

Communicating editor: S. L. Johnson 



604 I J. Liu eta/. 



:~-ZLG3' Genes | Genomes | Genetics 



